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I. INTRODUCTION 

B-meson decays to final states with baryons have been 
explored much less systematically than decays to meson- 
only final states. The first exclusively reconstructed de- 
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cay modes were the CLEO observations of B— >A+p7r 
and B— >A+p7T7r [l| and, later, of B°^D* + pp7r~ and 
B ^D*~pn 0. These measurements supported the pre- 
diction 3] that the final states with A c baryons are not 
the only sizable contributions to the baryonic i?-meson 
decay rate, and that the charm-meson modes of the form 
B—>D®NN'+anything, where the 2VW represent nucleon 
states, are also significant. Previous measurements show 
a trend that the branching fractions increase with the 
number of final-state particles. The branching fractions 
for the four-body modes B°^D < ^ + ppn~ @, Q are ap- 
proximately four times larger than those for the three- 
body modes B —^D^* >0 pp [51] , which, in turn, is approx- 
imately five times larger than those for the two-body 
modes B°^A+ p Q. 
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We expand the scope of baryonic _B-decay studies with 
measurements of the branching fractions and the kine- 
matic distributions of the following ten modes 0, [|[ : 



B°{ d 



Three-body B°-+D°pp 
Four-body B°^D + pp-K~ 

" B-^D°ppTT- 
Five-body B a ^D ppir^7r + 

" B-^D+ppn-ir- 



and B°^D*°pp, 
and B°^D*+ppir- , 
and B~^-D*°ppir~ , 
a,ndB°^D*°ppir-ir + , 
and B~^D* + ppTr~Tr~ 



Six of the modes — the four B~ and the two five-body B° 
modes — are observed for the first time. 

We reconstruct the modes through twenty-six decay 
chains consisting of all-hadronic final states (the list is 
given later with the results in Table [I]), e - g-j 

B + D*~ppTT + 7T + 
U D°7T- 

L> K + TT~TT + TT~ . 

A D° meson, as in the above example, is produced in 
eight of the B modes and a D + is produced in the remain- 
ing two. The £>°-meson candidates are reconstructed 
through decays to K~tt + , K~ir + ir°, and K ~ir + TT~ir + ; 
and the D+ to K~ "7r + 7r + . The D* -meson candidates 
are reconstructed through decays to D°ir° and the D* + 
asD°TT+. 

Typical quark-line diagrams for the three- and four- 
body modes with a meson are shown in Fig. [I] The 
three-body modes involve internal emissions of the W~ 
boson, whereas the four- and five-body modes involve 
internal and external emission diagrams. 

Baryonic B decays have a distinctive phenomenology 
whose features contrast with the patterns observed in 
meson-only final states. Experimentally, the overall rate 
enhancement of multi-body decays and the low-mass en- 
hancement in the baryon-antibaryon subsystem are ob- 
served [^4l^|. Theoretically, these modes are used to 
investigate a wide range of topics [l5- 25] . Among them 
are the predictions of the relative branching fractions, 
the decay dynamics, and the hypotheses involving exotic 
QCD phenomena, such as tetra-, penta-, or septa-quark 
bound states. In particular, there have been discussions 
of pp peaks near threshold values and penta-quark in- 
termediate resonance decays O c — >D® + p with respect to 
our modes [H-Sl]. 

The paper is organized as follows: Section HTl describes 
the data sample and the BABAR detector. Section IIIII 
presents the analysis method, introducing the key vari- 
ables Mes and AE. Section IIVI shows the fits to the 
joint Mes-AE 1 distributions. The fit yields and the cor- 
responding branching fractions are given. Section |V] dis- 
cusses the systematic uncertainties. Section IVII presents 
the kinematic distributions. For the three-body modes, 
the Dalitz plots of M 2 (D«°p) vs. M 2 (pp) are given as 
well as the invariant mass plots of the variables. For the 
four- and five-body modes, the two-body subsystem in- 
variant mass plots are given. In the four-body modes, we 
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FIG. 1: Typical quark-line diagrams representing (a) 
S°-KD wo pp and (b) B~—>D® pp7r~ modes. The gluon lines 
are omitted. 



investigate a narrow structure in the M(pir ) distribu- 
tion near 1.5GeV/c 2 . Section fVlII states the conclusions. 



II. BABAR DETECTOR AND DATA SAMPLE 

We use a data sample with integrated luminosity of 
414ft)" 1 (455xl0 6 BB) recorded at the center-of- mass 
energy ^=10.58 GeV with the BABAR detector at the 
PEP-II e + e~ collider. The e + and e~ beams circu- 
late in the storage rings at energies of 3.1 GeV and 
9 GeV, respectively. The value of ^/s corresponds 
to the T(4S') mass, maximizing the cross section for 
e+e~->- bb^T(AS)^BB events. The BB production ac- 
counts for approximately a quarter of the total hadronic 
cross section; the continuum processes e + e~— >uu, dd, 
ss, and cc constitute the rest. 

The main components of the BABAR detector [3(| 
are the tracking system, the Detector of Internally- 
Reflected Cherenkov radiation (DIRC), the electromag- 
netic calorimeter, and the instrumented flux return. 

The two-part charged particle tracking system mea- 
sures the momentum. The silicon vertex tracker, with 
five layers of double-sided silicon micro-strips, is closest 
to the interaction point. The tracker is followed by a wire 
drift chamber filled with a helium- isobutane (80:20) gas 
mixture, which was chosen to minimize multiple scatter- 
ing. The superconducting coil creates a 1.5 T solenoidal 
field. 

The DIRC measures the opening angle of the 
Cherenkov light cone, 9q , produced by a charged particle 
traversing one of the 144 radiator bars of fused silica. The 
light propagates in the bar by total internal reflection and 
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FIG. 2: Event display for the candidate decay B°-^D°pp, 
D°— >K + iv~ . The labeled tracks in the tracking system and 
DIRC rings at the perimeter correspond to the particles in the 
reconstructed decay chain. The remaining unlabeled tracks 
and rings are due to the decay of the other B° meson in the 
event. The beam axis is perpendicular to the image. 



is projected onto an array of photomultiplier tubes sur- 
rounding a water-filled box mounted at the back end of 
the tracking system. The DIRC's ability to distinguish 
pions, kaons, and protons complements the energy loss 
measurements, dE/dx, in the tracking volume. 

The calorimeter measures the energies and positions 
of electron-photon showers with an array of 6580 finely- 
segmented Tl-doped Csl crystals. 

The flux return is instrumented with a combination of 
resistive plate chambers and limited streamer tubes for 
the detection of muons and neutral hadrons. 

A data event display is given in Fig.[3]for the candidate 
decay B°^D°pp, W-^K+tt- . 



III. ANALYSIS METHOD 

This section describes the branching fraction measure- 
ment in four parts. Section IIII Al describes the Monte 
Carlo-simulated event samples that are used to evaluate 
the performance of the method. Section IIII Bl lists the 
discriminating variables and their requirements for the 
event selection. Section IIII CI defines the Mes and AE 
variables and presents their distributions for the newly 
observed modes. Lastly, Sec. IIII Dl describes the fit to the 
AIbs-AE distribution used to extract the signal yield. 



A. Monte Carlo-simulated event samples 

Monte Carlo (MC) event samples are produced and used 
to evaluate the analysis method. Two types of samples — 
signal and generic — are described below. 

The particle decays are generated using a combination 
of Evtgen [3l| and Jetset7.4 (32|. The interactions 
of the decay products traversing the detector are mod- 
eled by Geant 4 [33| . The simulation takes into account 
varying detector conditions and beam backgrounds dur- 
ing the data-taking periods. 

The signal MC sample is generated to characterize 
events with a B meson that decays to one of the signal 
modes (the accompanying B decays generically). The 
typical size of 3x 10 5 events per decay chain is two orders 
of magnitude larger than the expected signal in data. 

The generic MC sample is generated to characterize 
the entire data sample. The size is approximately twice 
that of the BABAR data sample. 

B. Event selection 

The e + e~ events are filtered for a signal S-meson candi- 
date through the pre- and the final selections. 

The pre-selection requires the presence of proton- 
antiproton pair and a D - or a _D + -meson candidate 
(written as D without a charge designation) in one of 
the 26 decay chains listed in Sec. Q] 

Protons are identified with a likelihood-based algo- 
rithm using the dE/dx and the 9q measurements as de- 
scribed in SecQl] For a 1.0 GeV/c proton in the lab frame 
(typical of those produced in a signal mode) , the selection 
efficiency is 98% and the kaon fake rate is 1%. 

The D-meson candidates are selected using the invari- 
ant mass [34[, M(D), and a kaon identification algorithm 
similar to that used for protons. The M(D) is required 
to be within seven times its resolution around the PDG 
value j.3j| (superseded later during final selection). For 
a 0.9 GeV/c kaon in the lab frame (typical of those pro- 
duced in a signal mode), the selection efficiency is 85% 
and the pion fake rate is 2%. 

For the D°^K-ir + n° and D*°^D°tt° sub-decay 
modes, the ir°— > 77 candidates are formed from two well- 
separated photons with 115<M(77)<150 MeV/c 2 or from 
two unseparated photons by using the second moment of 
the overlapping calorimeter energy deposits. 

The charged particles from the decay chain are re- 
quired to have a distance of closest approach to the beam 
spot of less than 1.5 cm. 

The final selection requires the presence of a fully- 
reconstructed signal B-meson candidate. Requirements 
on the discriminating variables described below are op- 
timized by maximizing the signal precision z—Sj V S+B 7 
where S is the expected signal yield using the signal MC 
sample and B the expected background yield using the 
generic MC sample. The signal is normalized using the 
measured branching fractions for the modes B a —>D®°pp 
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and B°—>D* + ppTr~ [1, @; for the rest of the modes the 
latter value is used. The quantity z is computed for 
each discriminating variable for each decay chain. For 
the variables with a broad maximum in z, the cut values 
are chosen to be consistent across similar modes. 

In order to select B-mcson candidates, M(D) is 
required to be within 3 (Jm(d) of the PDG value 
The resolutions a M{D) for D°^K-tt+, R-it+tt , 
K~tt + tt~tt + , and D + — >K~tt + h + arc approximately 6, 
10, 5, and 5MeV/c 2 , respectively. For the modes involv- 
ing D°— >K ~7T + 7r° decays, the combinatoric background 
events due to fake tt° candidates are suppressed using a 
model [36] that parameterizes the amplitude of the Dalitz 
plot distribution M 2 (K~n + ) vs. M 2 (7r + 7r ). The model 
accounts for the amplitudes and the interferences of de- 
cays of K*°^K-ir+, K*-^K-ir°, and p+^n+ir . The 
normalized magnitude of the decay amplitude is used to 
suppress the background events by requiring the quan- 
tity to be greater than a value ranging from 1% to 5%, 
depending on the mode. 

In order to select B*-meson candidates, the D*-D 
mass difference, AM=M(D°tt)-M(D°), is required to 
be within 3 a am of the PDG value [35|. The resolution 
a am is approximately 0.8MeV/c 2 for both D*°^D°ir Q 
and D*+^D°n + . For the mode B°^D°ppn-iT + , the 
requirement of AM>160 MeV/c 2 excludes the contami- 
nation from B°^D* + ppTT~ , D* + — >D°tt + decays. 

In order to select B-meson candidates, a combination 
of daughter particles in one of the signal modes is con- 
sidered. The momentum vectors of the decay products 
are fit [37j while constraining M(D) to the PDG value 
[35j . The vertex fit x 2 probability for non-B events peaks 
sharply at zero; these events are suppressed by requiring 
the probability to be greater than 0.1%. 

Continuum backgrounds events are suppressed by us- 
ing the angle ^thrust between the thrust axes [38[ of 
the particles from the B-meson candidate and from the 
rest of the event. The continuum event distribution of 
I cos ^thrust I peaks at unity while it is uniform for BB 
events, so the quantity is required to be less than a value 
ranging from 0.8 to 1, depending on the mode. 

After the selection, an average of 1.0 to 1.7 candidates 
per event remains for each decay chain and is largest for 
those decay chains with the largest particle multiplicity. 
If more than one candidate is present, we choose the one 
with the smallest value of 



(M(B) - M(B) PDG ) 2 (AM — AMpdg) 2 

5= ; - + ; - , (1) 



C. Definitions of M ES and AE 

The B meson beam-energy-substituted mass, Mes, and 
the difference between its energy and the beam energy, 
AE, are defined with the quantities in the lab frame: 



AE 




(2) 



The four-momentum vectors Qb—(Eb,Pb) and 
Qq=(Eo,Pq) represent the B-meson candidate and 
the e + e~ system, respectively. The two variables, 
when expressed in terms of center-of-mass quantities 
(denoted by asterisks), take the more familiar form, 

M ES =\/s/4- (Pb) 2 and AE=E*-^/2. 

The Mes-AE distributions for the events passing the 
final selection are given for the six newly observed modes 
in Fig. [3] Each point represents a candidate in an event. 
For many of the modes, a dense concentration of events is 
visible near Mes=5.28 GeV/c 2 , the PDG B-meson mass 
[HI, and AE=Q, as expected for signal events. The uni- 
form distribution of events over the entire plane away 
from the signal area is indicative of the general smooth- 
ness of the background event distribution. 

The M-es-AE plots are given in a box region of 
5.22<Af E s<5.30GeV/c 2 and |AB|<50MeV. This box 
is large enough to provide a sufficient sideband region 
for each variable where no signal events reside. It is 
also small enough to exclude possible contamination from 
other similarly related B-meson decay modes. 

For the purpose of plotting Mes and AE individually, 
the box region is divided into a signal and a sideband 
region. The Mes signal region is within 2.5<tm E s °f the 
mean value of the Gaussian function describing it and 
likewise for AE. Similarly, the Mes sideband region is 
outside 4 <7m es of the mean value and likewise for AE. 
The resolutions range from 2.2 to 2.5MeV/c 2 for ctm E s 
and 8 to lOMeV for <tae- The signal box is the intersec- 
tion of the Mes and the AE signal regions. 



D. Fit procedure 

The signal yield is obtained by fitting the joint Mes- 
AE distribution using a fit function in the framework of 
the extended maximum likelihood technique 39]. The 
likelihood value for N observed events, 



Wm{d)) 2 



(ctam) 2 



L(N,Cl) 



-N 



N 



Y[P( Vi ;N,n), 



(3) 



where the PDG values [35[ are labeled as such. The latter 
term in the sum is included only if a D* is present in the 
decay chain. If more than one candidate has the same S 
value, we choose one randomly. 



is a function of the yield estimate N and the set of pa- 
rameters f2. The yi is the pair of Mes and AE values 
for the B-meson candidate in the i th event and P is de- 
scribed below. The quantity L is maximized (40l - |42l ] with 
respect to its arguments. 



9 




FIG. 3: Scatter plots of Mes-A_E for the six newly observed B-meson decay modes: (a) B~~ — >D°ppn~ , (b) B°—>-D ppTT~iT + , 
(c) B~ — >D + pp-K~ n~ , (d) B~ — vD*°pp-K~ , (e) B°~^D*°pp-K~ii + , and (f) B~ ~^rD* + ppn~ -k~ . The first row of plots is related to 
the second by the exchange of the charm meson flf>D*. The decay chain involving D°— >K ~ty + or 7r + 7r + is shown. 

For (d, e, f), the decay chain involves D*°— >D°ir° or D* + — >D°7r + . The Mes projection, in lMeV/c 2 bins, is given above the 
scatter plot; the AE, in 1 MeV bins, on the right. For the projection plots, no selection is made on the complementary variable. 



The fit function is the sum of two terms 

P(yt\ N, Cl)=N sig P sig ( yi ; fl sig )+N hgd P hgd {yf, fi bgd ), (4) 

which correspond to the signal and the background com- 
ponent, respectively. For each component function, P a is 
the two-dimensional function, N a the yield, and Q a the 
parameters. The arguments of the function components 
are related to the quantities in Eq. ^ by N= J2p Np 

and f2= (J^ 

Each function component P a is written as the prod- 
uct of functions in Mes arL d AE since the variables are 
largely uncorrelated. (The signal bias due to the small 
correlation is treated as a systematic uncertainty.) The 
distributions for signal events peak in each variable, so 
P s i g is the product of functions composed of a Gaussian 
core and a power-law tail |43| . The background event 
distribution varies smoothly, so Pbgd is the product of 
a threshold function [38[ for Mes that vanishes at ap- 
proximately 5.29GeV/c 2 and a second-order Chebyshev 



polynomial for AE. 

The following function parameters are fixed to the val- 
ues found by fitting the signal MC distributions: the AE 
Gaussian width for P s i g , the Mes Gaussian width for Pig, 
the Mes power-law tail parameters for P s i g , and the A/es 
end-point parameter for Pb g d- Two exceptions are given 
after the detailed fit example. 

A detailed example of the fit results is given in Fig. [4] 
for the decay chain P~— >D + ppTr~TT~ , D + ^K~-k + tt + . 
The plots in Figs. 2^, andUJa show the Mes distributions 
for the AE signal and the AE sideband region, respec- 
tively. Likewise, Figs. St and |4jl show the respective 
AE distributions for the analogous Mes regions. The 
fit function projections describe the distributions in the 
sideband regions well (Figs. HJd, which gives us con- 
fidence that the background event distribution inside the 
signal box are also modeled well. 

The first exception to the fit procedure described above 
applies to the mode B~^D*°ppTr~ . A term is added to 
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FIG. 4: Fit details for B~ ->D + ppn~ n~ , D + ^-K~ty + ty + : 
Mbs and AE distributions in regions as noted on the plots. 
For (a, c) the top curve is the sum of P s i g and -Pbgd and the 
bottom curve is the latter; for (b, d) the curve is Pbgd- 



Eq. (j4]) to account for the sizable contamination from the 
mode B°—>D* + ppTr~. The fit function P poa k is the same 
form as P s i g with its parameters fixed to the values found 
by fitting the MC sample. The normalization iVp Ca k is 
based on the branching fraction measured in this paper. 

The second exception applies to four decay chains 
whose fits do not converge: P°— >D°ppTr~ir + , 
D°^K~ tt+tt ; B Q ^D ppir-ir + , D^^R-^tx-^; 
B ^D* ppn-7r+, P*°^P°7r n , D°^K-tt + tt°- 1 and 
B ^D* ppn-7r+, D*°^D 7r°, D ^K-tt+tt-tt+ . Two 
changes are made: the Gaussian parameters are fixed 
to the values found in the D°^K~tt + measurement, 
and the Mes end-point parameter is floated. The fits 
converge after the changes. 



IV. BRANCHING FRACTIONS 

This section presents the P-meson branching fractions B. 
Section llV Al shows the fits to the Mgg-AP distributions. 
Sections IIV Bl and IIV CI gives the B values and their ra- 
tios, respectively. Throughout this section, we simply 
state and use the systematic uncertainties of Sec. |V] 



A. Fits of M E s-AE distributions 

The Mes distributions for the events in the AE signal 
region for three-, four-, and five-body modes are given 
in Figs. [5HZ1 respectively. For all P-meson decay modes, 
the decay chains involving D ^K~it + or P + ^PT~7r + 7r + 
show a peak. 

The fit function projection in each plot describes the 
data well, except for the four decay chains corresponding 
to Figs. [7b, [7b, [7H, and[7f, which had difficulties with fit 
convergence as noted in the previous section. As we will 
see in Sec. [V] the yields from these decay chains do not 
contribute significantly to the P-meson branching frac- 
tion, which is dominated by the value from P°— >-PJ~7r + , 
because of their relatively large systematic uncertainties. 

The signal yields, given in Table HI range from 50 to 
3500 events per mode. 



B. Branching-fraction calculation 

The P-meson branching fraction for each of the twenty- 
six decay chains, given in Table U is given by 



B 



1 



1 



1 



2N B3 ByB D B D . e 



N sig - N, 



peak 



(5) 



whose ingredients are as follows: the number of PP pairs, 
7V B5 =455xl0 6 , the assumed T(45)^PP° or -^P+P~ 
branching fraction, By=l/2; the P-meson branching 
fraction. Bp the P*-meson branching fraction, 

Bd* [2a|; the reconstruction efficiency, e; the signal 
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FIG. 5: Mes fit projections for the three-body modes: (a-c) B°^D°pp and (d-f) B°^D*°pp, where (a, d) are reconstructed 
via D°^K-ir + , (b, e) D°^K--k+-k°, and (c, f) D° '-KfT ty+tt' 'tt 4 '; and (d-f) D*°^D°ty°. Events with AE within 2.5a of the 
mean value of the Gaussian function are shown. The top curve is the sum of P s i g and Pbgd and the bottom curve is the latter. 



yield, iV s i g ; and the measured contamination, A^ pea k, us- 
ing the M(D)-sideband data sample. The Bd' is in- 
cluded only when a D* decay is present in the decay 
chain. The efficiency e is determined using the signal 
MC sample and decreases with the particle multiplic- 
ity. The mode B°^D^pp, L»°^A'"7r+ has the highest 
value of e at 19% and B°~^D*°ppTr-n + , D*°^D°tt° and 
D°— >K~ir + Tr~ir + has the lowest at 1%. 

The B values, given in Table [Til are the combinations 
[HJ of the above measurements using the statistical and 
the systematic uncertainties. All B values are significant 
with respect to their uncertainties. For the previously 
observed modes, the results are consistent with earlier 
measurements . 



C. Branching- fraction ratios 



the charge of the charm mesons, D" + f>D^° , 

(iii) The ratio for the four-body mode to that of the 
corresponding three-body mode with one fewer pion is 
about four. 

(iv) The ratio for the five-body mode to that of the cor- 
responding four-body mode with one fewer pion is about 
one- half. 

The patterns (iii, iv) imply /?3-bod y <#5-bod y <#4-body 



V. SYSTEMATIC UNCERTAINTIES 

This section describes the systematic uncertainties for 
the B-meson branching fraction measurement. Sec- 
tion IV Al lists the sources, and Sec. IV Bl gives the error 
matrices. 



Table IIIII gives the ratio of the branching fractions B for 
modes related by D^D*, D® <&D® + , and the addition 
of n. These ratios show four patterns: 

(i) The ratios are roughly unity for the modes related 
by the spin of the charm mesons, D^D*. This result 
suggests that the additional degrees of freedom due to 
the D* polarization vector do not significantly modify 
the production rate. 

(ii) The ratio is roughly unity for the modes related by 



A. Sources 

The sources of systematic uncertainties, which are listed 
in Table ITVl can be organized as follows: 

• (i) Counting of the number of BB pairs, 

• (ii-iv) Assumed branching fractions, 

• (v-xi) Reconstruction efficiencies, 
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FIG. 6: Mes fit projections for the four-body modes: (a) B°^-D + ppTv , (b-d) B° — >D* + ppii , (e-g) B —>D°ppTT , and 
(h-j) B~^-D*°pp-K~ , where (a) is reconstructed via D + ->K~it + tv + , (b, e, h) D°-^K~-k + , (c, f, i) D° n + n° , and (d, g, 
j) D°-Kfr~7r + 7r _ 7r + ; and (b-d) D* + ^D°tx + and (h-j) D* ^D°tx . Events with AE within 2.5a of the Gaussian mean value 
are shown. For (a-g) the top curve is the sum of P s i g and Pbgd and the bottom curve is the latter; for (h-j) the middle curve 
is the sum of P pca k and Pb g d- 
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FIG. 7: Mes fit projections for the five-body modes: (a-c) B° — ¥D°ppn 7r + , (d-f) B°-^D*°pp-K ty + , (g) B —>-D + ppTT ir , and 
(h-j) B~ — >D* + ppir~ n~ , where (a, d, h) are reconstructed via D°— ¥K~n + , (b, e, i) D°— >K~ TT + ir a , (c, f, j) D° — >K~it + it~tt + , 
and (g) D + ^K~TT + ir + ; and (b-d) D* + ^D°ir + and (h-j) D*°->D 7r°. Events with AE within 2.5a of the Gaussian mean 
value are shown. The top curve is the sum of P s i g and the -Pbgd and the bottom curve is the latter. We note that the plots in 
(b, c, e, f) had difficulties with fit convergence; see text. 
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TABLE I: Intermediate values for TableHU B -meson branch- 
ing fractions for the decay chains. N B i s is the yield, -/Vpeak is 
the measured contamination (item xvii in Table ITV)) . and e is 
the reconstruction efficiency. The uncertainties are statistical. 
The rows marked by a dagger f have large systematic uncer- 
tainties; see text. The charges of the pions are implied as well 
as the D*°->D°7r and D* + ^D°n + decays, when applicable. 
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• (xii-xv) Fit functions and its parameters, and 

• (xvi-xvii) Backgrounds peaking in Mes or AE. 

These contributions are described below. 

(i) The number of BB pairs used in the analysis is the 
difference of the observed number of hadronic events and 
the expected contribution from continuum events. The 
latter is estimated using a separate data sample taken 
40MeV below the T(45) peak. The uncertainty of 1.1% 
is mostly due to the difference in the detection efficiencies 
for hadronic events in the data and the MC samples. 

(ii) The T(45) branching fraction is assumed to be 
equal for B°B° and B+B~. The uncertainty of 3.2% is 
the difference of 1/2 and the PDG value j3o| . 

(iii, iv) The D- and D*-meson branching fractions 
assume the PDG values [35[. The uncertainties of 
1.3%, 3.7%, 2.5%, and 2.3% are the PDG uncer- 
tainties for D°— S-if~7r + , K~n + ir°, K~ir + Tr~7r + , and 
D + ^K~n + n + , respectively; and 4.7% and 0.7% for 
D* a ^D°n° and D*+^D°tt + , respectively. 



(v) The charged track reconstruction efficiency is eval- 
uated using e + e~— >t + t~ events, where one tau decays 
leptonically and the other hadronically. The uncertainty 
of 0.5% is due to the difference between the detection 
efficiency in the data and the MC samples. 

(vi) The reconstruction efficiency of low-energy 
charged pion from D* + — >D°ir + decays is sufficiently dif- 
ficult, in comparison to other tracks, that item (v) cannot 
account for its uncertainty. Such a pion is often found 
using only the silicon vertex tracker because its momen- 
tum is relatively low. The momentum dependence of 
pion identification is evaluated using the helicity angle 
6?hei distribution — the angle between the pion direction 
in the D* + rest frame and the D* + boost direction — 
because the two quantities are highly correlated. Since 
the pions are produced symmetrically in cos#hci, the ob- 
served asymmetry in the distribution is indicative of the 
momentum dependence of the efficiency. The uncertainty 
of 3.1% is due to the difference in the momentum depen- 
dence in the data and the MC samples. 

(vii) The ir° reconstruction efficiency is evaluated using 
t + t~ events as in item (v) with an uncertainty of 3.0%. 

(viii) The signal B-candidate reconstruction efficiency 
is evaluated using the MC samples. Since these sam- 
ples use the uniform phase-space decay model while the 
reported baryonic decay dynamics 0, 0, [H, l9rjl4l this pa- 
per] are far from uniform, corrections are made in the 
variables where the strongest variation are seen — in bins 
of M 2 (pp) vs. M 2 (D^p)— using the data and the MC 
samples. The uncertainties ranging from 0.8% to 9.7% 
are due to the limited statistics of the samples. 

(ix) The particle identification efficiencies for kaons 
and protons are evaluated using the MC samples, which 
are then corrected using a data sample rich in these 
hadrons. The uncertainties ranging from 1.5% to 2.5% 
are due to the sample statistics associated with the cor- 
rection procedure. The sample, however, is dominated by 
the continuum events whose event topology is different 
from BB events. Items (x, xi) account for the differences. 

(x, xi) The kaon and proton identification efficiencies 
in the BB environment are evaluated using a data sample 
of D* + ^D°tt + , D°^-K~n + and A^pir~ decays, respec- 
tively. The uncertainties of 0.5% and 1.0%, respectively, 
are due to the differences in the event topologies. 

(xii) A subset of the fit function parameters is fixed 
when fitting the Mes~AE distributions in the data sam- 
ple. Such parameter values are obtained by fitting the 
MC distributions, and they are assigned an uncertainty 
from this fit. The effect on the signal yield is eval- 
uated by fitting the data sample with the parameter 
value shifted by lcr. The procedure is repeated for 
each parameter in the set. The uncertainties of 1.3%, 
2.8%, 5.7%, and 3.4% for the modes with D a ^R-n+, 
K~ir + n°, K~ir + n~-ir + , K~n + n + , and D + — >K ~ir + n + , 
respectively, are the quadrature sum of the fractional 
yield changes. 

(xiii) The choice of the signal fit function is evaluated 
using an alternate function, a fourth-order polynomial. 
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TABLE II: Main results of this paper: B-meson branching fractions for the ten modes. Also given are the values of , 
the degrees of freedom (DOF), and the \ 2 probabilities for the averaging of the results from Table [i] The measurements are 
consistent with the previous results. 



iV-body 


5-meson decay mode 


B ±0" s tat±0"syst 


X 2 /DOF Prob( X 2 ) 


R frnm Rpfs 9 ^ 


R frnm Rpf 4 












(96) 




ft r\~ 4 \ 

(1U ) 


Three-body B° - 


-> D°pp 


1.02±0.04 ±0.06 


4.3/2 


12 


1.18±0.15±0.16 [5] 


1.13±0.06±0.08 


a 


B°- 


D*°pp 


0.97±0.07±0.09 


4.1/2 


13 


1.20±g;ii ±0.21 [5] 


1.01±0.10±0.09 


Four- body 


B° - 


■> D + pp-iv~ 


3.32±0.10±0.29 








3.38±0.14±0.29 


// 


B° - 


•+ D* + ppiv~ 


4.55±0.16±0.39 


1.2/2 


54 


6.5 ±H ±1.0 [2] 


4.81±0.22±0.44 


// 


B~ - 


-¥ D°ppTT~ 


3.72±0.11 ±0.25 


3.4/2 


19 






// 


B~ - 


-)■ D*°ppn~ 


3.73±0.17±0.27 


0.5/2 


79 






Five-body 


B° - 


■> D°ppTT~ 7T + 


2.99±0.21 ±0.45 


0.3/2 


85 






// 


B° - 


■¥ D*°ppTT~n + 


1.91±0.36±0.29 


0.5/2 


78 






// 


B- - 


-¥ D + ppn~Tv~ 


1.66±0.13 ±0.27 












B~ - 


-> D* + ppTT~ TT~ 


1.86±0.16 ±0.19 


0.2/2 


91 







TABLE III: Ratios of B-meson branching fractions of the 
modes related by D ( * )0 <h>_D ( *) + , and the addition of 

7r. The uncertainties are statistical. 



Ratio of the modes 



R ± or 



Related 
B(B_° 
B(B° 
B(B~ 
B(B° 
B{B~ 

Related 
B(B~ 
B(ET 
B{B° 
B(B° 

Related 
B(B~ 
B(B~ 

Related 
B{B~ 
B(B~ 
B(B° 
B(B° 



by spin of charm meson 
-^D*°pp)/B{B Q ^_D°pp) 
->D* + ppn-)/B{B°^D + ppn-) 
-+D* a ppir )/B(B--*D ppir ) 
-^D*°pp-K--K + )/B{B ^D pp-K--K + ) 
- J >D* + ppn~n~)/B(B~ -^>D + ppn~n~) 
by charge of charm meson 
^D°ppn-)/B{B°^D + ppn-) 
->D*°ppir )/B(E°^D*' 
-^D°ppn-TY + )/B(B-^L 
-^D*°ppn-n + )/B(B~-^D* + ppTT-Ti 
by addition of pion to three-body modes 
->D*°ppn-)/B(B ^D* pp) 3 
->D°ppir-)/B(B°-*D pp) 3 
by addition of pion to four-body modes 
-^D + ppn'TT')/B{B°^D + ppn~) 
^D* + pp-K--K-)/B{B°^D* + pp-K-) 
-^D°ppn-n + )/B{B ^D ppn-) 



PpiT 



pp-K 



-■K + )/B{B°->D*°pp-K- 



95±0.08 
37±0.06 
00±0.05 
64±0.13 
12±0.13 

12±0.05 
82±0.05 
80±0.19 
03±0.21 

84±0.33 
64±0.18 

50±0.04 
41±0.04 
80±0.06 
51±0.10 



TABLE IV: Systematic uncertainty list for S-meson branch- 
fractions. The "D modes" represents D° — >K~tt + , 



mg 
K~ 



and K ty tt ty , and D + —tK 



Item Description 



Uncertainty (%) 



i Number of BB pairs 1.1 

ii R(T(4S): for T(4S)-kBB 3.2 

hi B(D): for D modes 1.8, 4.4, 3.2, 3.6 

iv B(D*): for D*-KD°7r°, D* + ^D°tt+ 4.7, 0.7 

v Charged particle reconstruction 0.5 

vi 7T+ from _D* + ^D°7r+ 3.1 

vii 7T° reconstruction 3.0 

viii Signal mode decay dynamics 0.8-9.7 

ix Kaon and proton id using data 1.5-2.5 

x Kaon id in BB event topology 0.5 

xi Proton id in BB event topology 1.0 

xii Fit function params: for D modes 1.3, 2.8, 5.7, 3.4 

xiii Signal fit function 0.6 

xiv Backgrnd. fit function: for D modes 0.8, 4.5, 1.3, 2.0 

xv Mes-A_E correlation 0.4-2.2 

xvi Backgrnd. peaking in Mes or AE for 0-5.5 (77-85) 
all modes (marked f in Table UJ) 

xvii Backgrnd. from baryonic modes 0.5-13.5 



The uncertainty of 0.6% is due to the yield difference 
with respect to the original fit function. 

(xiv) The choice of the background fit function 
is evaluated using a more general fit function with 
the addition of another such component. The 
uncertainties— 0.8%, 4.5%, 1.3%, and 2.0% for the 
modes with D°^K~ir + , K~ir + ir°, K~ir + ir~ir + , and 
D + ^K ~ir + n + , respectively — are due to the yield dif- 
ferences with respect to the original fit function. 

(xv) The small correlation between the Mes and AE 
distributions introduces a bias in the signal yield. This 
effect is quantified by fitting pseudo-experiments. Each 
experiment contains a background sample whose Mes 



and AE distributions are produced according to Pbgd, 
and a signal MC sample from the full detector simula- 
tion. The uncertainties ranging from 0.1% to 1.8% are 
from the deviation of iV s i g to the mean of the signal-yield 
distribution. 

(xvi) Background events whose distributions peak ei- 
ther at M E s=5.28GeV/c 2 or Ai?=0 can alter the signal 
yield. For the B~ ^D*°ppir~ measurement, the vari- 
ation of the normalization of the fit function for the 
B°^D* + ppir~ contribution within the experimental un- 
certainties has a negligible effect on the signal yield. For 
other B decay modes, no such sources are found. How- 
ever, the Mes distributions for a few cases feature a broad 
hump with a width around 20 MeV/c 2 spanning nearly 
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TABLE V: Systematic uncertainties (%) combined for the 
B modes. For each D mode, two columns are given. The 
uncorrelated values are given on the left columns and the 
correlated on the right columns. The right columns exclude 
items (iii, iv) of Table [TV] 



^^---^ — 


K 


_|_ 


K 7T 


-1- 


T7- — _| 


— 4 
TV 7T n 


K ir~ 


_ 4- 




unc 


cor 


unc 


cor 


unc 


cor 


unc 


cor 


B°^D°pp 


2.7 


3.5 


5.5 


6.9 


4.7 


7.0 






B°^D*°pp 


2.2 


4.6 


8.6 


8.6 


8.7 


7.7 






B°^D+ppn~ 














5.7 


5.7 


B°^D* + pp-K- 


4.2 


6.3 


7.6 


8.6 


9.9 


8.9 






B°^D°ppty-tt + 


14.5 


3.9 


81.2 


7.0 


77.3 


7.0 






B°^D*°ppn-n + 


13.8 


4.9 


86.3 


8.8 


85.4 


7.7 






B--^D°ppn' 


4.4 


4.2 


8.8 


7.2 


11.6 


7.3 






B-->D*°pp-n- 


6.6 


5.2 


8.4 


8.9 


10.5 


7.9 






B~^>-D + ppiT~n~ 














15.0 


5.9 


B~ — + ppTT~ 7T~ 


5.9 


6.8 


14.9 


9.0 


19.0 


9.2 







half of the signal box. The effect of the presence of such 
a source is quantified by adding a component P poa k to 
the fit function whose parameters are fixed except for 
the normalization. Except for four decay chains — those 
corresponding to Figs. [7b, Eh, Eh, and Ef — uncertainties 
ranging from zero to 5.5% are obtained from the changes 
in yield when the additional component is included. For 
the mentioned exceptions, the uncertainties range from 
77% to 85%. As a consequence of the large uncertain- 
ties, these four modes do not contribute significantly to 
the final results. 

(xvii) Background events from baryonic modes with- 
out a D meson are evaluated using the data sample. 
An example case where the final states are identical is 
B^Acpn , A c ^pK-n+ and B°^D°pp, D°^K-tt+t: . 
For such a source, the M(D) distribution does not peak 
at the D mass, so the contamination can be quantified by 
repeating the analysis with the M(D)-sideband region. 
^Vpeak is an additive correction factor for 7V s i g with un- 
certainties ranging from 0.5% to 13.5% due to the sample 
statistics. 



B. Error matrices 

The error matrix, V , spanning the D modes of a given 
B mode is the sum of the statistical and systematic com- 



■V, 



syst- 



ponents V= V si 

The Vgtat is diagonal with elements (cr s t a t,a) 2 (Table|I|. 
The V S y S t is the sum of a diagonal part and an off- 
diagonal part V Bya t=V }mc +V COI (Table E). The Kmc is 
diagonal with (a unCt0l ) 2 . The V cor is the sum of a diag- 
onal part with (<7 cor , a ) 2 and an off-diagonal part with 
Pa/3&cor,aCcor,i3- The correlation coefficient p a p is be- 
tween two D° modes a and j3. The correlations among 
D°-meson branching fractions are the PDG values [35| ; 
all others are assumed to be unity. 



VI. KINEMATIC DISTRIBUTIONS 

This section presents the kinematic distributions [45| . 
Sections IVTAl IVTBl and IVTCl give the plots for three-, 
four-, and five-body modes, respectively. Additional dis- 
cussion is devoted to the M(pw~) feature in Sec. IVIDl 

We briefly describe the background-subtraction and 
efficiency-correction methods used to obtain the differ- 
ential branching fraction plots (Figs. I MTTj) as a function 
of two-body invariant mass variables. The differential 
branching fraction, in bins j of the plotted variable, is 
the ratio of the number of signal events and the product 
of the correction factors as given in Eq. [5j The quantity 
in the numerator is the sum of the background-subtracted 
event weights for events in bin j; the formulae are given 
below. The efficiency-correction part of the denominator 
is found for bin j and is applied to each event weight. 

The S-Plot method is used [46| to find the event weight, 



5ig,sig-Psig(?/i) ~t~ Psig,bgd-Pbgd(yi 

N s - lg P sie {y t ) + N hgd P hgd ( yi ) 



(6) 



where the yi is the pair of A/es and AE values for the can- 
didate in the i th event; the fit functions P a were defined 
in Eq. (UJ). In general, the weight W is approximately 
for a background event and 1 for a signal event. The 
Psig.bgd quantifies the correlation between the signal and 
the background yields, 



N 

E 



Px{yi)Px>(.yi) 



Nsi g P s i g (yi)+N hgd P hgd (yi)Y 



(7) 



A. Three-body modes B^D^pp 

For the three-body modes, plots are given for Dalitz vari- 
ables and two-body invariant masses. 

The Dalitz plots of M 2 (D^°p) vs. M 2 (pp) for the 
events in the M^-AE signal box are given (Fig.[5h.,[Bb). 
The allowed kinematic region is the shaded contour. 

The background events present in Figs.[Shi and[U> are 
represented by Figs. |5t and [SU, respectively. The latter 
plots show the events in the MES-sideband regions with 
their normalizations determined from the background 
yield in the signal box. 

The two-body invariant mass plots are given in Fig. |21 
Differential branching fractions are plotted as a func- 
tion of M(D^°p) and M(pp) for events in different re- 
gions of the complementary variable. The two low-mass 
enhancements near threshold values in M(D"°p) and 
M(pp) correspond to the dense regions in the Dalitz 
plots. The broad enhancement in Fig. EJd) and |9jh) 
does not have a substantial contribution from J /tj) de- 
cays due to its width and current experimental limits on 
B°^D°J/i/}, J/i(>- 



47]. 



In general, we observe a strong similarity between the 
shapes of the corresponding distributions for B°^>-D pp 
and B°^D*°pp. 
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FIG. 8: Dalitz plots M 2 (pp) vs. M 2 (D"°p) for the three-body modes. Plots in the first column (a, c) correspond to B°— >D°pp; 
the second column (b, d) B°— >D*°pp. Plots in the first row (a, b) are the events in the Mes-AE signal box; the second row 
(c, d) the events in the MEs-sideband region normalized to the amount of background present in the respective plots in the 
first row. In the first row, near-threshold enhancements are seen compared to the respective sideband plots in the second row. 
The lines drawn at M 2 (pp)—5, M 2 (D°p)—9, and M 2 (D*°p) — 10.5 GeV 2 /c 4 are visual aides to show that the enhancements are 
mostly non-overlapping. The events are contained in the shaded contour representing the allowed kinematic region except for 
one outlier in (d), which failed the fit. The points are made larger for the plots in the second column for better visibility. 



B. Four-body modes B^-D^ppir 

For the four-body modes, plots are given for two-body in- 
variant masses in Fig. 1101 Differential branching fractions 
are plotted as a function of M(pp), M{D®p), M(D®p), 
and M(jpn~). 

The two-body invariant-mass distributions show a 



number of features. The M(pp) distributions show a 
threshold enhancement with respect to the expectations 
from the uniform phase-space decay model (Figs. [TUk . 
HDfc,II3,ITffln). The M(D®p) distributions show no in- 
dication of a penta-quark resonance at 3.1GeV/c 2 [H[ 
(Figs. [TOb. [Tot. [TO). [TOR The M(D®p) distribution in 
one of the modes (Fig.flOk) suggests a threshold enhance- 
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FIG. 9: Differential branching fraction plots for the three-body B-meson modes: (a-d) B° — >D°pf5 and (e-h) B°— ¥D*°pp. 
The captions give the various phase-space regions. The shaded region represents the uniform phase-space model with its area 
normalized to the data. The bin width for each row of plots is given on the left-most plot. 



merit, as was observed in the three-body modes, but the 
distributions in the other modes show no such features 
(Figs. [TOb, [TOfe, HOb). The M(pw~) distribution in one 
of the modes (Fig. UOt i) shows a narrow structure near 
1.5GeV/c 2 , but it is less prominent in the distributions 
of the other modes (Figs. [10h, US, HHp)- 

The peak near 1.5GeV/c 2 does not correspond to a 
known state. The peak is discussed in detail in Sec. IVI Dl 

C. Five-body modes B^D^ppirir 

For the five-body modes, plots are given for two-body in- 
variant masses in Fig. 1111 Branching fractions are plot- 
ted as a function of M(pp), M(D®p), M(D®p), and 
M(ptt"). 

In contrast to the distributions for the three- and four- 
body modes, the five-body distributions are generally 
more consistent with the expectations from the uniform 
phase-space decay model. 

A notable absence, again, is the signal of a penta-quark 
resonance at 3.1GeV/c 2 [48| (Figs. Eb, CUT, CCD, Eh)- 



D. Narrow M(pn ) peak at 1.5GeV/c 2 

The narrow peak in the M(jm~) [49( distribution at 
1.5GeV/c 2 , which we refer to as X, is discussed in this 
section. 

The opposite-sign M(p7r~) distributions corresponding 
to Figs. HUH . [TUh . ITU1 . [TUb are shown in more detail in 
Fig. HH In the detailed plots, the cc-axis bin width is 



smaller at 10 MeV/c 2 and the y-axis is the unweighted- 
uncorrected number of events. The events from the Mes~ 
sideband region is superimposed with its normalization 
determined from the background yield in the Mes-AP 
signal box. 

In order to measure the properties of the peak, the 
fit formalism of Eq. Q is used. The signal component 
P S i g is assumed to be a Breit-Wigner line shape. The 
background component -Pbgd is taken from the same- 
sign M(pn~) distribution. The distribution for the 
B°^D + ppir~ mode is relatively smooth (Fig. IT3a). and 
it describes the rise and fall of the opposite-sign distribu- 
tion well (Fig. IT^k ) , whereas the same-sign distributions 
in the other modes show a more rapidly falling behavior 
around 1.5 GeV/c 2 (Figs. H3b-[Hi). 

We note, however, that the use of the shape for Pbgd 
has limitations. Since the formation of the p or p is not 
necessarily symmetric with respect to the 7r~ in these de- 
cays, the same-sign M(pTr~) combination may not pre- 
dict the true shape for the non-resonant component in the 
opposite-sign M(pir~) distribution. As a consequence, 
we cannot precisely quantify the systematic uncertainty 
associated with the lack of knowledge of the true back- 
ground shape. 

For the two neutral B modes, the fits of the opposite- 
sign distributions describe the entire kinematic range well 
(Figs. H2k. [T2b). We note a small excess of events above 
1.65GeV/c 2 with respect to Pbgd, but no peak compo- 
nent is included in the fit at this mass. The fitted 
X mass is 1494.4±4.1 MeV/c 2 and 1500.8±4.4 MeV/c 2 , 
where the uncertainties are statistical, for B°— >D* + ppir~ 
and B a — >D + ppir~, respectively. We measure the full 
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FIG. 10: Differential branching fraction plots as functions of M(pp), M(D®p), M{D®p), and M(pn~) for the four-body B- 
meson modes: (a, b, c, d) B°^>-D + ppir~ , (e, f, g, h) B°— >D* + pp-K~ , (i, j, k, 1) B~ — >D ppn'~ , and (m, n, o, p) B~ ^D* pp-K~ , 
respectively. The shaded region represents the uniform phase-space model normalized to the data. The possible presence of 
a narrow peak near 1.5GeV/c 2 in plots in (d, h, j, p) are shown in detail in Figs. 112a . 112b . 112b . and 112b ?. respectively, and 
discussed in Sec. IVIDI The bin width for each row of plots is given on the left. 
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widths to be 51±18MeV and 43±17MeV, respectively. 
The widths are significantly wider than detector resolu- 
tion, which is less than 4MeV for a simulated X— >pn~ 
decay with a mass of 1.5GeV/c 2 and negligible width. 

In contrast to the neutral B modes, the opposite-sign 
distributions for the two charged B modes exhibit a less 
peaking behavior at I.5GeV/c 2 . As a result, the param- 
eter for the width in the B~^D°ppir~ mode is fixed to 
the value found in the B°^D + ppir~ mode; the results 
of this fit are not used in the average. 

The known nucleon resonances N* with the masses 
1440, 1520, 1535, and 1650MeV/c 2 are used in an at- 
tempt to describe the X. The distribution is fit with 
the N* fit function components each parameterized as 
a Breit-Wigner line shape. The normalization for each 
component is allowed to vary independently. However, 
the fit does not describe the peak because the X is much 
narrower than any of the N* resonances (Fig. IT4k). 

The overall significance of the X is difficult to mea- 
sure, due to our lack of knowledge of the true background 
shape, as discussed earlier, as well as further statistical 
issues. We caution that the X analysis is not blind, the 
parameters are not chosen a priori, and the distribution 
under the no-X hypothesis may be only approximately 
normal. Furthermore, even under the normal assump- 
tion, the presence of the mass and width nuisance pa- 
rameters under the alternative hypothesis means that the 
distributions of the S statistic is not likely to be pure x 2 - 

We provide a me asure of the statistical significance 
5=^/2 (In Li— In Lq) of the X in the two neutral B 
modes, where L\ is the likelihood value with P S j g and 
L is without P sig . The value is 5=8.6 for B°^D+ppir- 
and 5=6.9 for B°°^D*+ppn-. 

The systematic uncertainties are mainly due to the 
Pbgd- We fit using an alternate fit function by adding 
a component derived from the same-sign distribution of 
a different mode (Fig. 114b). The result is a mass shift of 
0.8 MeV/c 2 and a full width change of 4 MeV. An addi- 
tional contribution of 0.5 MeV/ c 2 is added for the mass 
measurement due to the absolute uncertainty of the mag- 
netic field and the amount of detector material [501 ] . 

In summary, the unknown structure X can be charac- 
terized by a Breit-Wigner line shape: 

M(X) = 1497.4 ± 3.0 ± 0.9 MeV/c 2 (8) 
T(X) = 47 ± 12 ± 4 MeV, 

where the uncertainties are statistical and systematic, re- 
spectively. 

VII. CONCLUSIONS 

We have presented a study of ten baryonic B-meson 
decay modes of the form B—tD^pptyfy using a data 
sample of 455 xlO 6 BB pairs. Significant signals are 
observed (Table HJ). Six of the modes — B~— >D°ppTr~, 
B-^D*°ppTT-, B ^D°ppTT-7r+, B°^D*°pptt-tt + , 



B -^D + ppir it , and B —>D* + ppTr tt — are observed 
for the first time (Figs. EHSfe, EMS), EH, EH, Ek, EM3, 
respectively). 

The B-meson branching fraction measurements range 
from 0.97xl0~ 4 to 4.55xl0" 4 with the hierarchy 
£>3-bod y <#5_body<#4-body ( Table HI)) . These results su- 
persede the previous BABAR publication of B Q ~^D°pp, 
D*°pp, D + ppn~ , and D* + ppn~ The branching frac- 
tions related by changes in the charge or the spin of the 
D meson are found to be similar (Table IIIII) . 

The kinematic distributions show a number of notable 
features. For the three-body modes, threshold enhance- 
ments are present in M(pp) and M(D®°p) (Figs, [g EJ. 
For the four-body modes, a threshold enhancement is ob- 
served in M(pp) and a narrow peak is seen in M(pir~) 
(Fig. HP]) . For the five-body modes, in contrast to the 
other modes, the distributions are similar to the ex- 
pectations from the uniform phase-space decay model 
(Fig. [TT]). 

The M(pn~) distributions in the neutral i3-meson de- 
cay mode show the most prominent peak near 1.5 GeV/c 2 . 
We obtained a mass of 1497.4±3.0±0.9 MeV/c 2 and a full 
width of 47±12±4 MeV, where the first uncertainties are 
statistical and the second are systematic, respectively 
(Figs. IT12HT4")) . Determining the significance and inter- 
preting the origin of the peak are complicated by the 
fact that the background fit function is parameterized 
by the distribution from the same-sign charge combina- 
tions pTT~ , a procedure which may not provide the true 
background shape. 

Despite the relatively small branching fractions for 
these modes of order 10 -4 , with product branching frac- 
tions of order 10~ 5 to 10~ 6 (including the D and D* 
modes), the large size of the BABAR data sample allowed 
us to observe signals containing hundreds of events in 
many of the modes. We are, therefore, able to probe their 
kinematic distributions that reflect the complex dynam- 
ics of the multi-body final states. 



VIII. ACKNOWLEDGMENTS 

We are grateful for the extraordinary contributions of our 
PEP-II colleagues in achieving the excellent luminosity 
and machine conditions that have made this work possi- 
ble. The success of this project also relies critically on the 
expertise and dedication of the computing organizations 
that support BABAR. The collaborating institutions wish 
to thank SLAC for its support and the kind hospitality 
extended to them. This work is supported by the US 
Department of Energy and National Science Foundation, 
the Natural Sciences and Engineering Research Council 
(Canada), the Commissariat a l'Energie Atomique and 
Institut National de Physique Nucleaire et de Physique 
des Particules (France), the Bundesministerium fur Bil- 
dung und Forschung and Deutsche Forschungsgemein- 
schaft (Germany), the Istituto Nazionale di Fisica Nu- 
cleare (Italy), the Foundation for Fundamental Research 



21 



(a) 



t )() 



2.0 2.5 3.0 

M(pp) (GeV/c 2 ) 




2.0 2.5 3.0 
M(pp) (GeV/c 2 ) 




2.0 2.5 3.0 

M(pp) (GeV/c 2 ) 




2.0 2.5 3.0 
M(pp) (GeV/c 2 ) 




3.0 3.5 4.0 
M(D°p) (GeV/c 2 ) 




3.0 3.5 4.0 

M(D*°p) (GeV/c 2 ) 




3.0 3.5 4.0 
M(D*p) (GeV/c 2 ) 




3.0 3.5 4.0 

M(D* + p) (GeV/c 2 ) 




3.0 3.5 4.0 
M(D°p) (GeV/c 2 ) 



1.5 2.0 

M(pjr) (GeV/c 2 ) 




3.0 3.5 4.0 
M(D*°p) (GeV/c 2 ) 



1.0 1.5 2.0 

M(pjr) (GeV/c 2 ) 




3.0 3.5 4.0 
M(D + p) (GeV/c 2 ) 



1.0 1.5 2.0 2.5 

M(pjr) (GeV/c 2 ) 




3.0 3.5 4.0 
M(D* + p) (GeV/c 2 ) 



1.0 1.5 2.0 

M(pjr) (GeV/c 2 ) 



FIG. 11: Differential branching fraction plots as functions of M(pp), M(D^p), M(D^p), and M(pn~) for five-body B- 
meson modes: (a, b, c, d) B°— >D°ppTT~ir + , (e, f, g, h) B°^>-D* ppir~TT + , (i, j, k, 1) B~ — >D + ppTT~ir~ , and (m, n, o, p) 
B~ — >D* + pp-K~-K~ , respectively. The shaded region represents the uniform phase-space model normalized to the data. Each 
B-meson candidate for the plots in (d, h) contributes two entries for both p-n~ and pir + combinations, so they are scaled 
accordingly. The bin width for each row of plots is given on the left. 
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FIG. 12: Fits of the opposite-sign M(pir~) distribution for (a) B°^D + ppir~ , (b) B° ^D* + ppir~ , (c) B~->D°ppir~ , and (d) 
B~ — >D*°ppTv~ for events in the signal box of Mes-A_E. The top curve is the sum of P a i g and -Pbgd while the bottom curve is 
-Pbgd- The Pbgd is from the corresponding plot in Fig. 1131 The shaded histograms are scaled Mes sidebands. A small inset plot 
is a close-up of the region around 1.5 GeV/c 2 ; its bin width is the same as in the larger plot. 
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